Valence bond crystal and possible orbital pinball liquid in a t2g model 
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We study a model for orbitally degenerate Mott insulators, where localized electrons possess t2g degrees of 
freedom coupled by several, competing, exchange mechanisms. We provide evidence for two distinct strongly 
fluctuating regimes, depending on whether superexchange or direct exchange mechanism predominates. In the 
superexchange-dominated regime, the ground state is dimerized, with nearest neighbor orbital singlets covering 
the lattice. By deriving an effective quantum dimer model and analyzing it numerically, we characterize this 
dimerized phase as a valence bond crystal stabilized by singlet resonances within a large unit cell. In the opposite 
regime, with predominant direct exchange, the combined analysis of the original model and another effective 
model adapted to the local constraints shows that subleading perturbations select a highly resonating ground 
state, with coexisting diagonal and off-diagonal long-range orbital orders. 
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I. INTRODUCTION 

One of the most intriguing concepts in contemporary con- 
densed matter physics is the possibility to stabiHze phases 
without any broken symmetry at low temperature. In strongly 
correlated materials, e.g. in transition metal oxides, charge 
degrees of freedom are localized due to large Coulomb inter- 
actions and both spin and orbital degrees of freedom interact 
with each other and may remain disordered. While numer- 
ous realizations of quantum spin liquids, such as resonating 
valence bond (RVB) phases, have been reported both experi- 
mentally and theoretically^^ there are only few examples of or- 
bital liquids in Mott insulators up to now.?— Usually orbitals 
are accompanied by spins and frustrated interactions in both 
subsystems could imply a spin-orbital liquid phase. If spin 
and orbital variables are coupled to each other by relativistic 
effects, such a phase may emerge from effective pseudospin 
interactions on a honeycomb latticei^ For the superexchange 
interactions in transition metal oxidesj^ a spin-orbital liquid 
could be stabilized in presence of geometrical frustration.- For 
instance, it was suggested as the ground state on a triangular 
lattice for LiNi02 with active eg orbitalsi^ It would explain the 
absence of symmetry breaking with spin and orbital order, but 
the realistic situation in this compound is more subtle]^ How- 
ever, it was shown recently that the paradigm of a spin-orbital 
liquid can be realized on the triangular lattice in a system with 
active t2g orbitals^li^ given by the (lll)-plane of Ti'^^ ions 
in NaTi02. 

The spin-orbital systems are rather complex and difficult 
to analyze as spin and orbital operators are frequently highly 
entangledii- and it is not trivial to recognize which part of the 
interactions is most frustrated. Therefore, considering sim- 
pler and easier tractable models with only orbital variables 



(and frozen spins) became recently fashionable as, in this 
way, the intrinsically frustrated orbital interactions could be 
directly investigated. These interactions have only the lattice 
symmetries,^ in contrast to high SU(2) symmetry of spin in- 
teractions. For a two-dimensional (2D) eg orbital model one 
finds competing yet robust orbital ordered phases jii In fact, in 
this case frustration is not yet maximal and increases further 
when the interactions are modified and the 2D compass model 
limit is approachedfi^ In this limit the system is still in the 2D 
Ising universality classic and orbital order persists in a range 
of finite temperature 

Three-dimensional (3D) orbital models are characterized 
by an even stronger frustration but an order by disorder mech- 
anism stabilizes an ordered phase in the Cg orbital model^^ 
as shown recently by Monte Carlo simulations 4^ The situa- 
tion in the 3D compass limit, also referred to as the classical 
t2g model, is still controversial: while the absence of a phase 
transition at finite temperature was conjectured by a high- 
temperature series expansionist a first order transition into a 
low-temperature lattice-nematic phase without any orbital or- 
der was found by Monte Carlo simulations^ 

Geometrically frustrated t-zg orbital systems on the trian- 
gular lattice have been studied recently by effective models 
for large Coulomb interactions."*"'^ These studies led to iden- 
tification of a variety of exotic ground states, depending on 
the parameters governing electron hoppings and interactions. 
Several of them were characterized by dimerization into ei- 
ther spin- or orbital-singlet phases — this enables a descrip- 
tion by an effective quantum dimer model (QDM) ;^-'^'^^ and 
opens a route towards the possibility of a spin-orbital liquidii 
in a regime with small Hund's exchange coupling where both 
spin- and orbital degrees of freedom can fluctuate. In a differ- 
ent situation where Hund's exchange coupling is large, a fer- 



romagnetic arrangement of spins is favored — here depending 
on which exchange mechanism is dominant, either an orbital 
singlet phase, or a non-dimerized phase characterized by lo- 
cal avoided-blocking constraints, were identifiedr^iii but sev- 
eral questions were left open about the precise nature of each 
phase. 

The aim of this article is to understand the effects of the 
competition between different exchange mechanisms natu- 
rally coexisting in such systems, by studying a microscopic 
orbital model and deriving effective Hamiltonians in extreme 
regimes. In Sec.|ll]we define the orbital model and analyze its 
low-energy spectrum and some observables to identify three 
regimes with distinct ground state properties. Then we ana- 
lyze the regime where superexchange dominates in Sec. HU] 
This leads us to derive and analyze, in SecHV] an effective 
QDM which allows us to characterize the corresponding phase 
as an orbital valence bond crystal (VBC). Later on, in Sec. |V] 
we focus on the opposite regime, where direct exchange dom- 
inates, and identify by several approaches an exotic phase, 
fluctuating and with long-range orbital order Eventually we 
summarize these findings in Sec. I VII 

II. MODEL AND GENERAL APPROACH 
A. Model definition 

We investigate an orbital model on a triangular lattice, with 
exactly one electron per site and t2g orbital degrees of free- 
dom. This model is a particular case of the one introduced in 
Ref. [13 to describe the limit of strong Coulomb interactions 
in layered, orbitally degenerate transition metal oxides, for in- 
stance the Mott insulating titanate NaTi02. In this compound, 
each titanium ion has a nominal valence 3+, corresponding to 
a single electron in the id shell; due to the crystal field result- 
ing from the octahedron formed by nearest neighbor oxygens, 
this electron is confined to the t2g subspace of this shell, char- 
acterized by the canonical basis: {\yz), \xz), \xy)}- In the 
present layered structure we can consider a single [111] plane 
of the NaTi02 crystal, with t2g electrons described by a model 
treating the combined effects of hoppings and Coulomb inter- 
actions on a triangular lattice. 

A general formulation of the electron dynamics and corre- 
lations requires a multiorbital Hubbard model,-^^!^ taking into 
account all t2g orbitals and spin degrees of freedom. For con- 
venience, in the following we label each bond direction with 
an index 7 e {a, 6, c} corresponding to the plane (respectively 
yz, xz and xy) in which the bond is embedded. Kinetic terms 
are of two types; (i) direct hopping between ions neighbor- 
ing on the triangular lattice, with amplitude t' — there, the 
electron can hop only between orbitals of a single flavor, de- 
pending on the direction of the bond [this flavor is e.g. xy = c 
for a bond in the xy plane of the cubic lattice, see Fig.[TJa,b)] 
— in consequence we rename the orbital flavor allowing for 
direct hopping in this direction with the same index,- that is 

\a) = \yz), \h) = \xz), \c) = \xy)- (1) 

(ii) indirect hopping with amplitude t, where an electron hops. 
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FIG. 1 : (Color online) (a) Representation of orbitals on two Ti ions, 
neighboring on a c-bond, and of the two distinct hopping processes 
between them: direct hopping of amplitude t' , and indirect hopping - 
via pz oxygen orbitals, with amplitude t. (b) An elementary triangle 
embedded in a [111] plane; each bond of this triangle has a label 
7 £ {a, b, c} which labels both the bond direction and [as indicated 
in Fig.[TJa)] the orbital flavor active in direct exchange on this bond. 



from a Ti ion, first to a neighboring oxygen and then to another 
Ti ion, neighbor of the former on the triangular lattice. This 
latter hopping process involves exclusively the two orbital fla- 
vors not involved in direct hopping on this bond; more pre- 
cisely, considering two ions neighboring on a c bond, an elec- 
tron initially in the a orbital of one ion can hop with indirect 
hopping only to the b orbital of the other ion, and vice-versa 
[see Fig.m^)]- Independently, the on-site Coulomb interac- 
tions are governed by Hubbard and Hund's exchange parame- 
ters U and Jh, respectively. 

In the limit where these interaction amplitudes are larger 
than hopping parameters, the id electrons localize and the 
system is Mott insulating. In this regime of parameters one 
can use an effective model acting within the subspace with 
one electron per site (in the present case of configuration 
of Ti'^+ ions), which describes interactions between their spin 
and orbital degrees of freedom. Such a model was derived 
for layered titanates in a previous work, see Ref. [l^, and con- 
tains exchange couplings. These follow from second-order 
processes where one electron hops from a site (Ti ion) to a 
neighboring site which becomes doubly occupied at an energy 
cost given by a linear combination of U and Jh, before one of 
the two electrons hops to the site left empty. The coexistence 
of direct- and indirect- hoppings implies that three exchange 
mechanisms are allowed. We list them below, indicating how 
their amplitudes depend on the hopping parameters t, t': 

(i) direct-exchange terms which result from processes where 
both hopping processes involved are direct hoppings, and are 
thus of amplitude cx t'"^ ; 

(ii) superexchange terms resulting from processes where both 
hoppings involved are indirect hoppings, and have an ampli- 
tude oc ; and finally 

(iii) mixed-exchange terms resulting from processes that com- 
bine one direct hopping and one indirect hopping, and their 
amplitude is thus oc ti! . 

These exchange terms couple both spin and orbital degrees of 
freedom of nearest neighbor (n.n.) ions on a bond {ij). As 
a consequence of the S'[/(2)-symmetry in spin space of the 
original Hubbard model, each of them can be written as the 
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sum of a term with a projection operator {Si ■ Sj — j) on a 
spin singlet state on the bond {ij), and a term with a projec- 
tion operator {Si ■ Sj + |), selecting a triplet state for this 
bond. 

We now specify further the context to a situation with large 
Hund's exchange coupling Jh (i.e., close to its maximal al- 
lowed value U/3), when only spin triplet states contribute to 
exchange processes. This favors the alignment of spins into 
a ferromagnetic phase, where exchange terms act only on or- 
bital degrees of freedom. This situation can actually also be 
realized, even if alone is not large enough to polarize 
spins, in the presence of an external magnetic field, the ef- 
fect of which adds to that of Hund's exchange — for NaTi02 
the estimated value of = Jh/U is ry ~ 0.14, i.e., not far 
from the value of the estimated transition to the spin-polarized 
phase rjc ~ 0.16;ii In such a spin-polarized phase, setting 
J= (i2+i'2)/([/_3Jj:^)astheunitofenergy(J= 1), the 
effective exchange Hamiltonian for orbital degrees of freedom 
is; 

n = {l- a)H, + ^/a{l - a)H,^ + aH^ , (2) 

where Hs, and i/m are, respectively, the superexchange, 
the direct exchange and the mixed exchange Hamiltonians de- 
fined below. The parameter 



interpolates continuously between the superexchange (iJg, for 
a — Q) and direct exchange (iJd, for a ~ 1) Hamiltonians — 
while the additional mixed-exchange term H-^^ is present only 
for < a < 1, in combination with the others two. The three 
terms of Eq. ^ are given by: 

= 2 ^ { {T+T+ + c.c.) + n,^n,^-,], (4) 
= - E E(^K7+H-c.) , (5) 

(»i>ll7 ^7^7 

i/d = 2 ^ Ui^Tij^. (6) 

{io)\\l 

The electron creation operator in orbital 7 = a, 6, c at site i 
is 7I, rii-f = 'yj^i, and the single-occupancy local constraint 
reads rii-y = 1. The notation (ij) || 7 indicates that a bond 
between sites i and j is oriented along the lattice direction 
7, see Fig. [Ub); in other words, this bond is parallel to the 
vector e-y defined by Cc ~ (1,0), e^/b = (±l/2,\/3/2) [see 
Figs. Oa) and [2b)]. In the expressions of super- and mixed 
exchange terms we use pseudospin operators Ti^ defined for 
each site and flavor index; for instance for 7 = 0; 

Tic = ^ ("m - nib) , (7) 

Ttc ^ 4K = {T^cY- (8) 

The operators T)a and Tib are obtained by cyclic permutation 
of flavor indices. The orbital flavor active in direct exchange 
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FIG. 2: (Color online) Clusters used for the numerical study of the 
orbital Hamiltonian Eq. (|2j, with; (a) iV = 16 sites and (b) iV = 12 
sites. The multiple occurrence of site indices indicates the periodic 
boundary conditions, i.e., invariance under translations by vectors 
Vi/2- (c) First Brillouin zone of the triangular lattice, with F, K = 
2Pa and Ma/b/c points shown — Mc = (0,27r/\/3) in terms of 
{kx,ky) coordinates. 

for each bond direction is depicted in Fig. [Itb); in superex- 
change on the same bond, both other flavors (and only those) 
are active. In Eqs. (jUi and (|5]l ^'j is the flavor index distinct 
from both /i and 7. 

B. Numerical analysis of the low-energy spectrum 

The Hamiltonian H [Eq. (|2]l] is studied by Lanczos exact di- 
agonalization (ED) on periodic clusters. Most generally, these 
are characterized by two (linearly independent) vectors Vi and 
V2; those with the smallest possible value of |Vi x V2I such 
that the image of any site by either of these vectors is identified 
to this site by periodicity. We choose to consider clusters in- 
variant under point group symmetries of the triangular lattice; 
for this, the simplest choice of vectors Vi is e.g. Vi = L.Cc, 
V2 = L.Ca, which coiTesponds to a = L^-site cluster as 
in the case N = 16, see Fig.|2ja). Another possible choice is 
to take Vi/2 = L.{ea + Cc/b)', the corresponding cluster has 
N = 3L^ sites; we will use hereafter the = 12 cluster 
shown in Fig.|2b) (and similar but larger clusters in Sec.lIVIi. 

The momentum-resolved low-energy spectrum gives the 
energy of the ground state, found at the F point. The a 
dependence of Eq, see Fig.|3la), shows a striking contrast be- 
tween the two opposite regimes; (i) a < 0.6(1) where Eq 
increases linearly with a, and (ii) a close to 1, where the a- 
dependence of Eq is oc Vl ~ Oi [see inset in Fig.Oa)]. The 
latter scaling indicates that the low-energy dynamics is here 
dominated by mixed-exchange terms. 

Considering low-energy excitations shown in Fig.|3jb), one 
notices an intermediate regime (iii) for 0.6 < a < 0.8, where 
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FIG. 3: (Color online) (a) Ground state energy per site Eq/N as 
function of a for periodic clusters of A'^ = 12 and 16 sites. Inset: 
square of the previous quantity as function of (l—a) in the vicinity of 
a = 1. (b) Lowest excitation energies AE{q) for momenta q = F, 
Ma, and (depending on cluster size) either K or (tt, 0), see Fig.[2tc). 



several quasi-degenerate lowest states at T and points are 
well separated from higher excitations. We checked that there 
are exactly six such low-energy states. Three of them, at the F 
point, consist of a pair of exactly degenerate states, which are 
the exact ground states for a finite range of a in this regime,— 
and a third low-energy state which is the lowest otherwise. 
The boundaries ctcr.i of the zone with exact ground state de- 
generacy depend slightly on cluster size, e.g. the upper bound- 
ary is acr,+ ^ 0.81(1) and 0.72(1) for iV = 12 and iV = 16 
respectively. The three remaining low-energy states are de- 
generate, one being at each M-^ point. These features indi- 
cate a gapped, ordered phase with six-fold degeneracy of the 
ground state in the thermodynamic limit (TL). 

In contrast, in regime (ii), although lowest states are also 
found at F and points indicating a symmetry-broken 
phase, above them, excited states of energy ~ c^^g^l — a are 
compatible with modes becoming gapless in the TL — for de- 
tails see discussion in Sec.[V] Eventually, for < a < 0.6(1), 
the lowest excitations are at the F point and have energies 
0(1) decreasing with increasing N, compatible with an order 
breaking non-translational symmetries. Note, however, that 
since both clusters considered are of moderate size and with 
different shapes, we cannot characterize this phase from their 
low-energy spectra alone. 



III. DIMERIZATION IN THE SUPEREXCHANGE 
REGIME 

In this Section we focus on the regime where superex- 
change is the dominant exchange mechanism (i.e., we con- 
sider the regime of a ^ 1). We will provide, within the orbital 
model Eq. (|2|i, evidence for a dimerized ground state, which 
we will study further in the next Section using an effective 



FIG. 4: (Color online) (a) Dimer occupation number nd{c), and (b) 
the resonance amplitude P4 [Eq. (TTJ, with site indices as indicated 
on the left-hand-side lozenge; each lozenge in the inset shows a sin- 
glet covering at play in the resonance] for A*" = 12 and 16 as func- 
tions of a. Shaded areas indicate the range of possible values due to 
an exact twofold ground state degeneracy occurring for both clusters 
over a finite range of a - see text and Ref.lTTI. 



model. 

The superexchange interactions alone, in the present con- 
text, are known to favor dimerization into orbital singlets on 
nearest neighbor bonds. This is clear when one considers 
an isolated bond {ij) \\ 7, for which the ground state of 
is a singlet wave function in terms of the pseudospin variables 
Ti^^ and Tj,~^. On larger systems, we will evaluate the strength 
of dimerization by computing the dimer or singlet expectation 
value, defined as; 



where I^Po) is the ground state found in ED, and 
1 



b]h\ 



(9) 



(10) 



is the operator creating a singlet from the electron vacuum for 
a c-bond {ij) \\ c. While for a ^ 1 one finds rather small val- 
ues of Ud — 1/9, see Fig.|4ja), a much larger value ~ 0.36(2) 
is found for a <C 1, indicating a dimerized ground state. Note 
that the value Ud — 1/9 is close to the values expected for 
uncorrelated orbitals and almost unaffected by the avoided- 
blocking constraints. Note as well that the shaded areas in 
Fig. E] for 0.6 < a ^ 0.8 are not resulting from computa- 
tional uncertainty, but from the exact ground state degeneracy 
occurring in the corresponding (size-dependent) range of pa- 
rameters. Depending on the actual ground state |4'o) selected 
within the ground state manifold (here a 2D space), for an ob- 
servable O not commuting with T-L, (5'o|C l^'o) can take all 
possible values in the range indicated by the shaded area (and 
delimited by the same symbols as in the non-degenerate case). 

Interestingly, the value of nd obtained for a <C 1 is close 
to 7/24 obtained for a variational state I'l'var), built as an 



equal-amplitude superposition of the twelve columnar singlet 
coverings which minimize the energy among 
static singlet coverings. However, the true ground state has a 
much lower energy and singlet correlations significantly less 
modulated than those estimated with l^'var), see Appendix A. 
Hence we compute the quantity: 



P4 = 



{d{^dl^di2d34 + H.c.) 



(11) 



defined on a lozenge and shown as function of a, see Fig.|4jb). 
A discussion of the value of P4 on an isolated lozenge and its 
interpretation in terms of singlet resonance are given in Ap- 
pendix A. The large values of P4 shown in Fig. Hfb) for the 
dimerized ground state of the superexchange regime suggest 
that it is stabilized by such resonances between nearest neigh- 
bor singlets. These resonances can favor either a valence bond 
crystal (VBC) or a dimer liquid; in the first case large quan- 
tum fluctuations may reduce strongly the amplitude of any or- 
der parameter related to this order, in comparison with the 
amplitude expected for a model VBC wave function (such as 
I'I'var) for a columnar phase). Thus, the identification of the 
phase cannot be addressed directly in ED within the orbital 
model Eq. (|2]l, but requires to use an adapted effective model, 
which we describe hereafter 



IV. QUANTUM DIMER MODEL 
A. Derivation of a Quantum Dimer Model 

Motivated by the evidence discussed above for a dimer- 
ized phase in the superexchange regime, we derive a quantum 
dimer model (QDM) to provide a better understanding of this 
phase. Such a model, thanks to the reduced number of de- 
grees of freedom at equal size (here the degrees of freedom 
are not anymore orbital configurations at each site, but the po- 
sitions of dimers in close-packed dimer coverings fulfilling a 
hard-core constraint), can be accessed numerically for much 
larger system sizes than the original orbital model of Eq. (|2|. 
It is then possible to address the issue of the behavior in the 
TL. For this purpose, we follow the Rokhsar-Kivelson (RK) 
schemer^is^ in the present context this means that we project 
the orbital Hamiltonian in the subspace spanned by nearest 
neighbor orbital singlet coverings \C). An important charac- 
teristic of these coverings is that they form an overcomplete, 
thus non-orthogonal basis. Thus, instead of diagonalizing the 
Hamiltonian projected onto the singlet coverings' subspace 
(with matrix elements {C'\T-l\C)) one has to solve a gener- 
alized eigenvalue problem involving both this matrix and an 
overlap matrix with elements (C \C) . Due to the high connec- 
tivity of these matrices the numerical solution of this problem 
is, a priori, much more time consuming than the diagonaliza- 
tion of a sparse matrix of the same size. 

An approximation done in similar cases22 consists of ex- 
panding the matrix elements mentioned above in powers of 
a parameter x - here, this parameter is defined formally 
such that the singlet wave function on a bond (ij) \\ c is 
x{\aiaj) — \bibj)) (similar expressions hold for singlets on 
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FIG. 5: (Color online) Kinetic processes (three left columns) and 
configurations with finite potential energies Ac — ±v (two right 
columns) in the effective QDM Hamiltonian Hqdm of Eq. il2i de- 
rived in the superexchange regime. The coupling between configura- 
tions in the leftmost column, with amplitude t4, is found at 4*'' order 
in overlap expansion, with x = l/\/2; subsequent columns show 
couplings on two inequivalent six-bond loops, with amplitudes tea 
and tai, found at sixth order in this expansion. A potential term with 
amplitudes ±v applies to every interdimer bond when dimers touch- 
ing this bond are parallel to each other (the sign depends on whether 
they are parallel or not to the interdimer bond). 



bonds oriented along either Ca or Ch). Note that, although x 
takes the value 1 / in the case of interest, here it is intro- 
duced as a control parameter for the perturbative expansion. 
An effective Hamiltonian matrix, in an orthogonal basis of 
dimer coverings | c) which are in one-to-one correspondence 
with singlet coverings |C), is then obtained at given order p of 
this expansion. 

We carried out this derivation up to order p = 6, and ob- 
tained the following effective Hamiltonian: 

■Hqdm = ^ A,|c)(c| -J2tc,c'\c){c'\ , (12) 



where dimer coverings |c) and |c') differ by a shifted closed 
loop of length I = 4 or 6. The off-diagonal terms of ampli- 
tudes tec' and diagonal terms of amplitudes are explicited 
below and schematized in Fig. |5j these amplitudes are ex- 
pressed in units of the superexchange amplitude J(l — a). 

(a) Off-diagonal (kinetic) terms shown in three first 
columns of Fig. |5] — each term cx |c)(c'| couples all 
pairs of dimer configurations which differ from each 
other only by the position of dimers on a short loop 
of length I < 6. These terms originate from the off- 
diagonal part of superexchange interactions on inter- 
dimer bonds. The term with amplitude ti ~ — 4a;^, 
which flips two dimers on a ^ = 4-long loop, is found at 
fourth order in the derivation; at sixth order, two terms 
flipping three dimers on loops of length / = 6 have am- 
plitudes tga = 3t4/4 and igf, = ^4/4, depending on the 
shape of the loop. 

(b) Diagonal terms cx |c)(c|, which play a role of effec- 
tive potential terms (found at order 4), act on interdimer 
bonds and favor certain energetically optimal states. 
These terms originate from the diagonal part of su- 
perexchange terms on these bonds. The resulting po- 
tential energy Ac can be written (after a global energy 
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FIG. 6: (Color online) Real-space dimer correlations obtained for the 
QDM Eq. ( 112b on periodic clusters: (a,b) on the N = 36 cluster, for 
(a) w — 1.5 and (h) v = 0.5; and (c) on the = 48 cluster for 
V = 0.5. In (a) correlations correspond to the ground state of the 
(1,1) topological sector only, see Sec. lIVC] 



shift) as a sum over interdimer bonds of terms taking 
values ±v when both dimers touching a bond are par- 
allel to each other (see Fig. |5]l, and otherwise. In the 
present derivation we have found that v = 1^/2. 

This effective QDM mimics the accurate properties of the 
original model in this regime, similarly as in Refs.iandi^for 
spin-singlet phases. Note that the potential energy is different 
from the one of the RK QDM on the triangular lattice24 (where 
it is proportional to the number of lozenges with dimers oc- 
cupying two parallel edges); here the potential acts on each 
interdimer bond instead of lozenges. 

Although we derived Hqdm for the purpose of studying the 
superexchange limit a = 0, it is also relevant in the regime of 
< a ^ 1 where dominant superexchange couplings are per- 
turbed by small mixed-exchange couplings (with amplitude 
~ y/a). Indeed, matrix elements of the latter do not con- 
nect different singlet coverings, thus they do not modify the 
effective QDM (the direct-exchange perturbation bring some 
potential terms which, if taken into account, would modify 
^QDM, but should be of amplitude cx a and thus negligible in 
the limit considered). The main effects of mixed- and direct- 
exchange terms are to contribute to quantum fluctuations out 
of the n.n. singlet covering subspace, and to destabilize the 
dimerized phase for increasing a. 



B. Analysis of the quantum dimer model 

We have studied Hqdm using two numerical methods: (i) 
ED on periodic clusters with = 12, 36 and 48 sites — 
for the largest one the use of translational symmetries and 
topological invariants (see Sec. llVCl l allows us to reduce the 
Hilbert space size to Nh ~ 1.0 x 10^ in a representation of 
momentum F; and (ii) a zero temperature Green's function 
quantum Monte Carlo (GFMC).^*' This method can be per- 
formed here since all the off-diagonal terms are negative. It 
allows us to obtain the ground states of significantly larger 
clusters than within ED; here we focus on periodic clusters of 




FIG. 7: (Color online) (a) Dimer-dimer order parameter M(q) 
[Eq. ( 1131 )1 q = Pa or q — Ma [see Fig.|2lc)] as function of u ob- 
tained from ED on clusters of A = 12, 36 and 48 sites, (b) Scaling 
of M{q) (obtained from GFMC) with inverse linear size for v — 0, 
and either q — Ma (circles) or q — Pa (squares); dashed lines are a 
guideline to the eye. 



A^ = 3L^ sites up to L = 10. We choose to set, unless they 
are explicitly specified, the amplitudes of off-diagonal terms 
to their values obtained in the derivation, and adopt t4 as unit 
of energy in this section. But we keep ?; as a free parameter — 
this choice allows us to compare the case of interest {v = 0.5, 
corresponding to the superexchange limit a = of the orbital 
model) to two limiting cases that are easier to characterize, 
and which we address first: (i) v ^ 1 and (ii) \v\ <C 1. 

In case (i), for = the potential energy is mini- 
mized by 0{2^) degenerate ground states (maximally fiip- 
pable states2i). Once quantum fluctuations are turned on 
via tc.c' 7^ 0, this degeneracy is lifted and a particular 
ordered phase is selected by a quantum order-by-disorder 
mechanismi^^ for \teb\ < l^eal one finds^i a columnar VBC. 
In particular, for igf, = iea/S = 1/4, this is confirmed by 
the real-space correlations depicted in Fig. |6ja) on a 36-site 
cluster and at u = 1.5; these correlations clearly indicate a 
columnar pattern. 

More quantitatively, the associated translational symmetry 
breaking is reflected by a Bragg peak, at points q = of the 
Brillouin zone, in the dimer-dimer order parameter 



M{q) = ^ { ( *o 



E 



1/2 



$0 



C13) 



where |$o) is the ground state for a cluster considered. Here 
dj,e;, ~ (1), if a dimer is absent (present) on the bond 
(ik) II 7 such that the image of the site i by the translation 
of is k. For v > 1.5 and A^ > 36, M{Ma) exceeds 95% 
of the maximal value 1 / a/T2 obtained for a fluctuation-free 
columnar order, see Fig. Ha); and this order parameter stays 
unambiguously finite in the TL for v > 1.0(2). 

For case (ii), dimer correlations suggest a different transla- 
tional symmetry breaking in the TL, shown in Fig.|6jb,c). The 
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Bragg peaks in M (q) found at both Ma and Pa [see Fig. 13 a), 
and size scaling for v = in Fig.|7tb)] clearly indicate an or- 
dered pattern invariant by 27r/3-rotations, with a 12-site unit 
cell called a plaquette in the following. This phase has been 
found to be stabilized by an Z = 4-loop kinetic term, either 
isolated^i or along with moderate I = 6-loop terms bring- 
ing extra resonances within a plaquettei^ When v increases, 
AI{Pa) decreases gradually, in parallel with the increase of 
M{Ma), and eventually vanishes in the TL in the range where 
the large values of M{Ma) indicate a columnar order, see 
Fig. |2ta)- For a value v > 0.5, one observes in this figure 
pronounced jumps in both order parameters, and between this 
value and v ~ 0.8 it is not clear whether both quantities, or 
only M{Pa), vanish in the TL. This can correspond either to 
the onset of a RVB spin liquictfi or to a transition point be- 
tween two distinct VBCs. Unfortunately, for v > 0.6 the 
GFMC suffers from a lack of convergence, which restricts the 
available cluster sizes,-2iS2, while the correlations obtained by 
ED are more consistent with a columnar phase. 




FIG. 8: (Color online) Topological gap AEt [Eq. l il4b 1 versus v ob- 
tained for the QDM using large clusters of A'^ = 3L^ sites with peri- 
odic boundary conditions. The data are obtained by exact diagonal- 
ization for L = 4 and with quantum Monte Carlo for L > 6. 



C. Topological gap 

We also discuss the behavior as function of v of another 
quantity easily obtained within the effective QDM: the so- 
called topological gap A£'t^22i2^i2^ which is well defined a on 
periodic system, e.g. a torus, for a QDM where off-diagonal 
terms correspond only to local updates (this is the case here). 
On a torus accommodating the rotation symmetries of the lat- 
tice, it is defined as follows, 

AEt = \E„{Wc = Wa^l)- EoiWc = Wa = -1)1 , 

(14) 

where = (—1)"^ are topological invariants, such that 
is the number of dimers crossing a line parallel to the 7 axis 
and winding around the torus. Since the parity of Ua and ric is 
conserved by all kinetic terms of the QDM, each topological 
sector, or subspace spanned by all dimer configurations with 
fixed Wa and Wc, is well-defined. We also remark that, on 
clusters invariant by point group symmetries of the triangular 
lattice, the two sectors such that Wa + Wc = have the same 
spectrum as either the Wa = Wc = — 1 sector, or the Wa = 
Wc = 1 one, depending on the cluster size. 

The topological gap provides valuable insights into the na- 
ture of the ground state. In a VBC this gap grows as the linear 
size L of the system. This is because the ordered pattern fits in 
only one of the two topological sectors {Wa , Wc) — (1,1) and 
(—1,-1); and AEt corresponds to an excitation disturbing 
the crystalline order over a whole winding loop. On the con- 
trary, in a dimer liquid all topological sectors are fully equiv- 
alent in the TL and the topological gap is simply a finite-size 
effect, e.g. for a gapped Z2 liquid it behaves as L^e~^/^, with 
^ > and 5 a constant depending on the parity of L/2.— 

In the present case, on a large cluster of size N ~ [see 
Fig. [8] with L < 8], AEt decreases linearly with v from a fi- 
nite value at u = [expected in a plaquette phase where the 
ground state is in the sector (Wc = Wa = (—1)'^^^)] down 
to almost zero for Vc — 0.55(5). When v increases further 



it takes again a finite value in the columnar phase: here, the 
ground state is found in the sector {Wc = Wa = 1) (and, 
if L/2 is odd, also in sectors such that Wc = —Wa). Thus, 
when going from one phase to another, depending on the par- 
ity of L/2 the ground state changes or not its topological sec- 
tor; we have verified that the value of v minimizing AEt {v) is 
roughly independent of L, so that this minimum is a good indi- 
cator for the transition. In contrast, the maximum of AEt{v) 
for V ~ 0.9 in the = 48 case merely signals a crossing 
between a non-local excitation of energy AEt oc (as 
expected from perturbation theory in the large v regime) and 
another non-local excitation which is of lower energy close 
to the transition to the plaquette phase, where this excitation 
energy vanishes. From this data the hypothesis of an interme- 
diate RVB phase, in which At extrapolated to the TL would 
vanish in an finite range of v, seems unlikely. Besides, the 
vanishing of the topological gap observed for Vc — 0.55(5) 
coincides with the changes occurring in the dimer order pa- 
rameter. 

Although we cannot exclude a liquid phase stabilized in a 
narrow range of v ~ 0.5, the most probable scenario is thus a 
first order transition between plaquette and columnar phases, 
occuring at w ~ 0.55(5); this implies that the former phase is 
expected to be stabilized in the superexchange regime. 



V. THE DIRECT EXCHANGE REGIME 

We concentrate now on the regime a > 1 /2 of the model, 
where direct exchange is dominant. For a ^ 1 all other cou- 
plings can be treated as perturbations and, as we did in the 
previous sections, it is again possible to derive an effective 
Hamiltonian which correctly captures the ground state prop- 
erties of the microscopic system. In the following, we report 
the presence of an exotic ground state in this regime, and pro- 
pose a complete analysis for its characterization. 



8 




FIG. 9: (Color online) Orbital states for a ~ 1: (a) typical mixed- 
exchange process connecting two avoided-blocking states, which dif- 
fer on the shaded bond; (b) ordered patterns — the purely coUinear 
pattern (top) is frozen for a — 1, while the other (bottom) allows 
for low-energy dynamics on the highlighted bonds by consecutive 
processes triggered by mixed-exchange terms (arrows), whereas one 
sublattice (marked by full dots) is ordered. 



A. Effective model within the avoided-blocking subspace 

We start from the limit a = 1, in which H = Hd is 
positively definite and selects a macroscopically degenerate 
ground state manifold, characterized by rii^rij^ = on every 
bond {ij) parallel to e^. These orbital configurations, called 
avoided-blocking states, can be described in a representation 
where a rectangle at each site of the triangular lattice repre- 
sents the orientation of the occupied orbital at this site — a 
rectangle parallel to 7-bonds (but centered on the site i, in 
contrast to dimers of Sec. |IV]which were centered on bonds) 
means that Ui^^ = 1 in the configuration considered here, see 
Figs. |9ja) and|9lb). The number of such configurations grows 
exponentially with system size, i.e., they are macroscopically 
degenerate and form a low-energy subspace separated from 
higher-energy configurations by a gap ~ 2q! ~ 2 in this limit. 
Close to this limit (for 0<1 — a^l) the weight of other 
configurations, with blocking defects — consisting on bonds 
{ij) II 7 where on both sites i and j the 7 orbital is occupied 
— is proportional to 1 — a. This scaling can be easily under- 
stood within second order perturbation theory where mixed- 
exchange terms, which can create/annihilate such defects, are 
treated as perturbations w.r.t. direct exchange. 

In this regime, the ground state is thus mostly composed 
of avoided-blocking configurations, and its properties are de- 
termined by the action of mixed-exchange (and to a lesser 
extent of superexchange) perturbations within the avoided- 
blocking subspace. Indeed, as shown in Fig. |9|a) a single 
mixed-exchange term can connect distinct avoided-blocking 
states, allowing for important quantum fluctuations. We in- 
vestigate whether fluctuations can account for a disordered 
ground state, or whether the perturbing exchange terms select 
a long-range orbital order, by a quantum order-by-disorder 
mechanism. To that extent, we define and consider an effec- 
tive Hamiltonian, found at first order in a perturbation theory 
where super- and mixed exchange are considered as perturba- 
tions to the direct exchange Hamiltonian: 

H^s = VHsV + VUn^V) . (15) 



Here V is the projection operator onto the avoided-blocking 
subspace, i.e., in the basis of orbital configurations {|c)} 
(eigenstates of all "fjj.i operators), Vic) = |c) if |c) is an 
avoided-blocking configuration, while 7^|c) = otherwise. 

The ground state |^o*) of Hcs can be accessed for signifi- 
cantly larger cluster sizes (N < 36 sites in Lanczos ED) than 
the ground state \'^q) of the original Hamiltonian H; thus it 
allows us to address in a more controlled way the possible 
existence of long-range order in the TL. 

B. Orbital ordering 

The first quantity we consider to address the possibility of 
orbital ordering is the structure factor of orbital correlations, 

^->^^ -^T. ^'■^'-'^'^ \n^,n,,\ VI/) , (16) 

choosing for \ '^) either the ground state |4'o) of the full model 
or an approximate state I^Pq*) of the effective model defined 
in Sec. IV Al Due to local constraints imposed by dominating 
direct-exchange interactions, Sc{q) is expected to be maxi- 
mal at Ala and Mb. The type of ordered pattern which max- 
imizes this quantity consists of lines with orbital flavors uni- 
form within a line and alternating between neighboring lines 
— we call this collinear order. On a cluster invariant by point 
group symmetries, Sc{Ma) = 1/12 if 1^1*) in Eq. (fT6] l is cho- 
sen as an equal weight superposition of all such patterns (there 
are six of them, since lines parallel to e-y and containing 7 or- 
bitals are forbidden here). Here we find that for a larger than 
0.6(1), Sc{Ma) takes values of the same order of magnitude, 
and the comparison between different clusters [see Fig. fTOl' a)] 
indicates that in the TL this quantity may stay finite. 

In the limit a 1 on which we focus first, this is con- 
firmed by Sc{Ma) obtained from the ground state I^Pq*) of 
the effective model Heft, which for a — > 1 reduces to — 
Hcs{a ^ 1) cx VHraV; Fig. [Totb) clearly indicates that 
Sc{Ma) computed with I^Pg*) stays finite in the TL, scaling 
roughly as ci + C2/N where ci > 0. Note that its value for 
e.g. iV = 12 matches well the one in Fig. fTOl a) for the same 
cluster size and a — > 1; besides, both wave functions \'^o) and 
l^tg*^) have there a large mutual overlap ~ 1 — ka. This con- 
firms the validity of the effective model Eq. ( fTSb in this limit, 
where the ground state is thus characterized by long-range or- 
bital order inducing a translational symmetry breaking in the 
TL. 

A collinear phase would be compatible with such a re- 
sult: indeed the corresponding ordered patterns [see Fig.|9tb)- 
top] belong to the avoided-blocking subspace. Yet, the stabi- 
lization of this phase is not expected in this limit; collinear 
ordered patterns do not allow for mixed-exchange fluctua- 
tions, whereas these are essential in selecting a ground state 
with finite and negative (iJm) (within the effective model we 
estimated {H^)/N ~ -0.35(5) in the TL). With a simi- 
lar reasoning, one can exclude other candidate phases with 
static orbital order. A dimerized phase could also be con- 
sidered as a candidate, since mixed-exchange favors a res- 
onating state on individual bonds, e.g. the resonating state 
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>^N=16 
■N=12 




O-O^l/N 0-08 



FIG. 10: (Color online) Structure factor Sc(Ma) of orbital correla- 
tions, see Eq.lll6t. as obtained (a) within the orbital model of Eq. (|2} 
for the ground states of AT = 8, 12, 16 clusters, and for increasing a; 
and (b) within the effective model of Eq. dlSb for the ground states of 
N < 36 clusters, for a — > 1; for computational details see Ref . [sol. 
In (b) the dashed line is a guide to the eyes. 



(1 / \/2){\ciaj) + \biCj)) on a bond (ij) \\ c; yet such a phase, 
where blocking defects cannot be avoided, would not be fa- 
vorable either. 

While the nature of the ground state for a ^ 1 remains 
open after the discussion above, one notices that an upturn 
occurs in Sc{Ma) when a is decreased until reaching a value 
ac ~ 0.81(1) (ac ^ 0.86(2)) for iV = 12 (A^ = 16), see 
Fig. [iota). This upturn coincides (N — 12) or not (TV = 16) 
with the boundary of the regime with exact twofold ground 
state degeneracy; it signals a change in the ground state, which 
for a < ac could be collinear-ordered. Indeed, the diagonal 
part of superexchange is minimized among avoided-blocking 
states by collinear patterns, in which e.g. a and b orbitals 
are never neighboring on a c-bond; far enough from the a ^■ 
1 limit, superexchange becomes nonnegligible compared to 
mixed exchange and can favor the collinear order. 



C. Insights from low-energy excitations 

We now analyze, still in the direct-exchange-dominated 
regime, the low-energy spectra of both the orbital Hamilto- 
nian H [Figs. [TTT a) and [TTT b)! and the effective Hamiltonian 
Hctf [Figs. [TTT c) and [TTT d)1. For a close to 1, as mentioned 
in Sec. IIIBI the lowest states are found in representations of 
momenta F and AI^.— Above them, one notices excitation 
branches with energy 



(17) 



i.e., proportional to the amplitude of mixed-exchange pro- 
cesses. In contrast to low-energy states characteristic of a hy- 
pothetical gapped, orbital ordered phase evoked earlier, these 
states are not restricted to high-symmetry momenta, and thus 
they (or at least some of them) are not a mere consequence 
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FIG. 1 1 : (Color online) Top part — low-energy spectrum of the or- 
bital Hamiltonian "H in the direct-exchange-dominated regime (a > 
0.8) for clusters of (a) iV = 12 and (b) A*" = 16 sites. Bottom pai't 
— low-energy spectrum of the effective Hamiltonian Hch [Eq. l ll5b 1. 
for a > 0.7, and clusters of: (c) iV = 24, and (d) iV = 16 sites. 



of translation symmetry breaking in the TL. When comparing 
between different sizes the values minq-{c7v.,j'} coiTesponding 
to the lowest of these branches (which is typically at F and 
points) this value decreases significantly with increasing 
and seems to vanish in the TL; excitations with a similar 
scaling law [Eq. ( fT7] )l are also found at other momenta — they 
have energies getting smaller with increasing A^. These states 
may follow a certain dispersion law Ld{q), quantized here by 
the finiteness of clusters considered, but corresponding to a 
gapless mode of this phase. 

This situation contrasts with what we see further away from 
the direct-exchange limit, that is, for 0.6 ^ a < 0.8 within 
the effective model — as well as within the full model for the 
same range of a, see Fig.[3lb). In the latter regime, the spec- 
trum is characteristic of the collinear phase with lowest states 
only at F and Mj, while excitations with other momenta are 
found at much higher energies. Between both regimes, a level 
crossing occurs between two groups of quasi-degenerate low- 
est states, of different nature — on each side of the crossing 
the number and momenta of lowest states indicate the sym- 
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FIG. 12: (Color online) (a) Amplitude of on-site orbital fluctuations 
(fofc) and (b) von Neumann entanglement entropy Svn{A), both 
evaluated in the ground state of Eq. ^foi N — 12, 16 site clusters^, 
as a function of a. The inset shows the subset A of lattice sites used 
in the definition of Svn{A) [see Eq. il9\]. 



metries broken in the TL. At equal size, this crossing occurs 
for a value of a slightly larger in the original model than in 
the effective model (e.g. for N = 16, at a ~ 0.9 and ~ 0.8, 
respectively). The origin of this difference in the crossing po- 
sition is that the effective model underestimates effects of off- 
diagonal superexchange terms which contribute to stabilize 
the collinear phase. But the crossing position is not much sen- 
sitive to the system size (within the effective model, on clus- 
ters = 12, 16, 20, 24, this crossing occurs at a = 0.755(5), 
0.795(5), 0.73(1), and 0.79(1) respectively). This shows that 
the effective model Eq. (fTsT i captures well, qualitatively, the 
transition between the collinear phase found in a regime where 
super-, direct and mixed exchange amplitudes are compara- 
ble, and a phase with a distinct, unconventional orbital order, 
found in the direct exchange regime i.e. for a close to 1 . 



D. Strong orbital fluctuations in a symmetry-broken phase 

The unconventional structure of the low-energy spectrum 
observed for < 1 — a < 0.2, with possibly gapless exci- 
tations, contrasts with the — at least partial — orbital order- 
ing evidenced by large values of Sc{Ma)- To understand how 
these features can coexist, we evaluate, on the ground state of 
the full model Eq. (|2]i and for a fixed site i, the amplitude of 
on-site orbital fluctuations: 



(6tc) = ( 



(18) 



shown in Fig. [T2? a). This quantity is remarkably large for 
a > 0.9, where it has only weak size-dependence. We 
checked that on-site orbital fluctuations have comparable am- 
plitudes for larger clusters, if computed within the effective 
model Eq. ( fTSl l in the same range of a. The importance of 
these fluctuations is consistent with a finite orbital compress- 
ibility, i.e., the fact that the electron density in 7 orbitals can 



be tuned continuously by an adapted chemical potential (for 
details see Appendix B); and with the presumably gapless ex- 
citation modes discussed in the previous paragraph. 

When a is decreased further away from the direct-exchange 
limit, a drop in (6^c) occurs simultaneously with the upturn 
in Sc{J\Ia) and the crossing in the low-energy spectrum dis- 
cussed previously. This is a further indication for the transi- 
tion toward static collinear order for a smaller than this value, 
while for a close to 1 the mixed-exchange energy stabilizes 
a qualitatively more fluctuating ground state. From the data 
of Figs. |3] [To] and[TT] one can locate the transition between 
these two phases at roughly a ~ 0.8(1). 

Another instructive quantity is the von Neumann entan- 
glement entropy between one sublattice and the rest of the 
system,— defined as: 



Svn{A) = -Tr {pA In(pA)} , 



(19) 



where pA is the reduced density matrix computed for the 
ground state of an A^-site cluster^^ corresponding to a par- 
tition between an ensemble A of N/ 4 sites forming a triangu- 
lar sublattice with doubled unit cell [see inset of Fig. fT2l b) or 
Fig.|9jb)-bottom] and the complementary sites, which form a 
kagome lattice. This quantity, shown in Fig. [TSlb). displays 
a clear jump at the same position (on the a axis) as the jump 
in (fe^c) discussed previously. The entropy iS„jv(^) is signif- 
icantly larger on the right hand side of the jump, evidencing 
here a more fluctuating ground state, than on the left hand 
side. Besides, we compared this quantity with entanglement 
entropies corresponding to other partitions; in those cases, we 
also found jumps at the same value of a, but with smaller 
amplitudes [e.g. when considering Svn{B) with B being a 
lozenge of 4 sites, on the A^ = 12 cluster, between a = 0.8 
and a = 0.82 this quantity increases by ~ 0.15 compared to 
an increase of ~ 0.47 in Si,n{A)]. 

These features, in a ground state with spontaneously bro- 
ken translation symmetry, motivate us to suggest the follow- 
ing scenario: the symmetry breaking corresponding to trans- 
lations by vectors (7 = a, b, c) allows the ground state 
to have a structure where four different sublattices play non- 
equivalent roles; in a simplified picture [see Fig.|9tb)-bottom] 
one sublattice is ferro-orbitally ordered, with e.g. c-orbitals 
covering this sublattice. On other sites, which alone form 
a kagome lattice, the number of electrons of a given orbital 
flavor fluctuates mainly via mixed-exchange processes. By 
these processes, an electron with c orbital flavor can propa- 
gate along c bonds as long as it doesn't meet another elec- 
tron of this flavor; at each step of this propagation, the elec- 
tron moving in the opposite direction converts its orbital fla- 
vor from a to 6 or vice versa. In other words, such a succes- 
sion of mixed-exchange processes on neighboring bonds (see 
plain and dashed arrows in Fig. |9lb)-bottom) can be seen as 
the effective motion, at arbitrary distance and without direct- 
exchange-induced energy cost, of a single c orbital along a 
line containing alternating a and h orbitals. Similar effective 
single-orbital motions are also possible in other directions; the 
low-energy orbital dynamics are thus not confined in one spa- 
tial dimension but rather on an effective kagome lattice. The 
freezing of the complementary sublattice may be interpreted 
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as an order-by-disorder effect. One can also note that the lat- 
tice symmetry breaking evidenced here has been predicted for 
interacting t2g electrons in such a geometry, due to the struc- 
ture of hopping s. -"* 

Altogether, this scenario allows to understand qualitatively 
that, in the direct-exchange-dominated regime for a > 0.8, 
the ground state is characterized by: (i) large on-site orbital 
fluctuations, (ii) a significantly higher entanglement entropy 
S-uN (A) than in the collinear state favored for smaller a, and 
(iii) above the former ground state there may be a continuum 
of low-energy excitations, possibly gapless. This phase would 
be an orbital analog of a (bosonic) supersolid^l or (fermionic) 
pinball liquid,^ and is in any case an original and rather exotic 
type of ground state in the context of orbital models. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have addressed a situation where large 
Coulomb interactions lead, even though the spin degrees of 
freedom are frozen in a ferromagnetic phase, to a high frus- 
tration of effective interactions between t2g orbital degrees of 
freedom. This is due to two factors: first, even when con- 
sidered separately, both superexchange or direct exchange in- 
teractions are frustrated on the triangular lattice, while ana- 
log models on e.g. a square lattice would allow for orbitally 
ordered ground states; second, in the present model there is 
a competition between three distinct exchange mechanisms 
which naturally coexist in the compounds such as NaTi02, 
motivating this model. To address the nature of the ground 
states selected by these interactions, we studied both the or- 
bital model itself and two effective models adapted to extreme 
situations. For these studies we employed numerical tech- 
niques — mostly exact diagonalization, extended by Monte 
Carlo techniques for one of these models. This allowed us to 
characterize several ground states, depending on the param- 
eter a governing the ratios between amplitudes of the three 
exchange interactions: 

(i) In the regime of dominant superexchange {a <C 1), we 
found concordancing evidence for a dimerized phase, where 
electrons pair into orbital singlets on nearest neighbor bonds. 
The effective quantum dimer model which we derived in this 
context allowed us to characterize this phase as a plaquette va- 
lence bond crystal (VBC) with a large unit cell (of 12 sites), 
and which is stabilized by resonances between singlets within 
this unit cell. We also found that slight modifications of inter- 
action parameters might lead to a transition to another, colum- 
nar, VBC; this suggests that such a system should remain dis- 
ordered down to very low temperatures. Yet, the tendency 
towards dimerization and the mechanism stabilizing the pla- 
quette VBC seem to be robust to small mixed exchange inter- 
actions perturbing the superexchange ones, for a < 0.6(1). 

(ii) In the intermediate regime where superexchange, mixed 
exchange and direct exchange interactions have comparable 
amplitudes, for 0.6(1) ^ a < 0.8(1), we identified a differ- 
ent phase, with collinear orbital ordering in lines, where the 
orbital flavor is uniform within a line but alternates between 
neighboring lines. This phase, with gapped excitations, seems 



to be stabilized by the joined effects of direct exchange and 
superexchange interactions — the former select a large num- 
ber of low-energy orbital configurations {avoided-blocking 
states), among which those favored by superexchange cou- 
plings are the collinear patterns. 

(iii) Eventually, in the third regime, for a > 0.8(1) where 
direct exchange is by far the dominant exchange mechanism, 
we found an original type of ground state, resulting from the 
action of subdominant mixed exchange interactions within 
avoided-blocking states. This ground state is characterized by 
a spontaneous symmetry breaking, which allows for a struc- 
ture where orbital long-range order and strong orbital fluc- 
tuations develop on distinct sublattices; the corresponding 
phase has low-energy modes, presumably gapless, within the 
avoided-blocking subspace. The spatial ground state struc- 
ture indicated by our results suggests that it might be an origi- 
nal realization of supersolid- or pinball-like behavior in orbital 
physics. 

Based on these findings, several open questions are left for 
future studies: in the superexchange-dominated regime, the 
proximity to the columnar-plaquette transition within the ef- 
fective dimer model we considered indicates that a refined 
dimer model, including terms found at further order in the 
overlap expansion, might predict a columnar ground state, or 
even a third dimerized phase, for the superexchange regime. 
In the direct-exchange-dominated regime, it would also be 
valuable to confirm — or invalidate — our scenario by com- 
plementary techniques, ranging from mean field approaches 
to variational techniques, applied either to the orbital model 
or to the effective model adapted to this regime. Indepen- 
dently from these issues, one could also investigate whether 
the phases we evidence here might be realized in existing 
compounds — for this purpose one should include small octa- 
hedral distortions (which exist in NaTi02 iIjtIi ) along with the 
various exchange couplings considered here. An investigation 
of properties of the resulting orbital model, if compared to hy- 
pothetical (e.g. thermodynamic, or of resonant inelastic x-ray 
scattering) measurements on titanates in high magnetic fields, 
would allow to estimate the relative amplitudes of the different 
exchange mechanisms; one could then consider these ampli- 
tudes (or this value of a) and focus on the low-field regime, to 
study the possibility of a spin-orbital liquid phase. 
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Appendix A: Resonance processes in the superexchange limit 

In the superexchange regime, for < a <C 1, the tendency 
toward dimerization is quite strong and favors a priori a large 
number of nearest neighbor orbital singlet coverings, where a 
singlet is created (from the vacuum) by the operator dlj on a 
bond (ij) — for instance, for a bond (ij) \\ a, 



(Al) 



Distinct singlet coverings arise naturally due to the form of 
Hs and are non-orthogonal to one another. In a dimerized 
phase, the ground state composed of these coverings may un- 
dergo resonances which lower its energy and determine the 
orbital correlations. To quantify the importance of resonance 
processes between different singlet coverings we define the 
quantity P4 related to resonances on the smallest non-trivial 
loop, that is on a lozenge. By labeling sites on a lozenge in a 
way shown in Fig. |4] this quantity reads as: 

Pi = (^O d\^dl^di2d34: + 4244^^13^24 ^o) , (A2) 

where |5'o) is the ground state found in exact diagonaliza- 
tion on the cluster considered. Typical values of P4 are 
P4 ~ 0.29(2) in the superexchange limit a = 0, as seen in 
Fig. lUb), and decrease with increasing a in the same way 
as the dimer expectation value rid- This indicates that, while 
fluctuations out of the dimerized subspace increase, the reso- 
nance processes stabilizing the dimerized phase subsist until 
this phase becomes less favorable than a non-dimerized phase. 
The values found on = 12, 16 clusters can be compared to 
the extreme (maximal) value corresponding to a 2 x 2 clus- 
ter, where such a resonance (or flip) affects the whole system. 
Assuming here again periodic boundary conditions (so that 
all orbital flavors are equivalent), one can easily check that 
the ground state of this cluster at a = reads: 

|^'o"°) = ■^||aia2a3a4) + I&162&3&4) + |ciC2C3C4)| 

^ --j=^\aia2b3b4) + \b1b2a3a4) + \a1C2C3a4) 

+ \c1a2a3C4) + \b1C2b3C4) + \cib2C3b4)Y (A3) 



and from this it follows that 



P4 = - 



V2 
3 



0.971 



(A4) 



Such a large value cannot be expected on significantly larger 
clusters, where resonances on lozenges coexist with reso- 
nances on loops of length / > 6; but even there, values 
P4 ~ 0.29(2) for a — > prove that the contribution of such 
resonances to orbital dynamics is significant in this regime, 
leading to the ground state energy per site Eq ~ —1.25(5). 
Note that this energy is much lower than expected for a static 
singlet covering. In fact, the superposition l^'var) of columnar 
singlet coverings has an energy per site. 



Appendix B: Influence of an external chemical potential in the 
direct-exchange regime 

An interesting aspect shedding light onto ground state prop- 
erties for a — > 1 is its sensitivity to an external orbital field, 
or chemical potential, distinguishing one type of orbital (say, 
c ones) from the two others. The effective model augmented 
by this potential reads now: 



cff.c 



off + M< 



(Bl) 



Similarly, within the full orbital model, one can add such 
a chemical potential of amplitude pc- When this one is 
positive and large compared to the mixed-exchange ampli- 
tude Xm = y/a — a^, the chemical potential selects, among 
avoided-blocking configurations, those which minimize (to 
zero) the population of c orbitals defined as 



(B2) 



-0.5. 



(A5) 



The selected configurations are two collinear states with al- 
ternating rows of a and b orbitals, see Fig. |9lb)-top. For 
jic ^ —Xjn (but both quantities small before the direct- 
exchange amplitude), the chemical potential favors the four 
other collinear states, obtained from the two previous ones by 
a lattice rotation combined with a flavor permutation, and now 
maximizing Uc within the avoided-blocking subspace. We 
show in Fig. [13] the dependence of Uc on /^tc within the con- 
strained model Hcs.c in the direct-exchange limit a 1, as 
well as within the fuU orbital model for two values, small but 
finite, of 1 — a. 

Instead of a direct transition between the two ordered 
regimes mentioned above, characterized respectively by Uc = 
and 71c = 1/2 up to quantum fluctuations, we observe an 
intermediate regime where ric varies continuously with Hc- In 
this regime, the ground state has a finite orbital compressibil- 
ity (defined here as ^^)- When the mixed-exchange ampli- 
tude Xm is small, the values of the chemical potential corre- 
sponding to transitions from the intermediate regime to both 
collinear ordered ones, iic,± [see Fig.[T3]for /^c,+], are close 
to ±Xm- When a is decreased, one sees that pc.+ decreases 
slightly — this behavior is consistent with our indications for 
a collinear phase when /ic = and a < 0.8. Results for the 
constrained model confirm the presence of the intermediate 
regime, with a finite negative slope of nc{pc) in the vicinity of 
Pc = 0, separated from a laige-pc regime where ric = by a 
succession of discontinuities in the range 1 < /ic/ x'm < 2 (the 
multiplicity of discontinuities is a finite-size effect). The finite 
orbital compressibility at and in the vicinity of pc = 0, com- 
bined with the long-ranged orbital correlations [evidenced e.g. 
by Sc{Ma) in Fig. (TO) constitute similarities between this 
phase and a bosonic supersolid: in both situations, order de- 
velops on one sublattice, while either the number of bosons or 
the population of c-orbitals can fluctuate on the other sublat- 
tices, allowing for a finite compressibility. 
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FIG. 13: (Color online) Electron density in c orbitals per site, ric, 
as a function of the ratio of over mixed-exchange amplitude 
Xm = ^/(l — ct)a, either within the full model {H augmented 
by the chemical potential term) on periodic clusters (A^ = 12, 16 
sites) for different a values; or within the constrained model H^s^c 
[Eq. l IBlH . as obtained for clusters of A'^ < 24 sites for a — 1. 
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